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ABSTRACT 

We investigate the effect of modified gravity on the specific angular momen- 
tum of galactic halos by analyzing the halo catalogs at z = from high-resolution 
A^-body simulations for a f{R) gravity model that meets the solar-system con- 
straint. It is shown that the galactic halos in the f{R) gravity model tend to 
acquire significantly higher specific angular momentum than those in the stan- 
dard ACDM model. The largest difference in the specific angular momentum 
distribution between these two models occurs for the case of the isolated galac- 
tic halos with mass less than 10^^ h~^MQ, which are likely least shielded by the 
chameleon screening mechanism. As the specific angular momentum of galactic 
halos is rather insensitive to the other cosmological parameters, it can in principle 
be an independent discriminator of modified gravity. We speculate a possibility 
of using the relative abundance of the low surface brightness galaxies (LSBGs) as 
a test of GR given that the formation of the LSBGs befalls in the fast spinning 
dark halos. 

Subject headings: cosmology:theory — methods:statistical — large-scale struc- 
ture of universe 

1. INTRODUCTION 

It has been almost one century since Einstein introduced his theory of general relativity 
(GR) as a mathematical framework within which the evolution of the Universe can be coher- 
ently described. Remarkably successful as it has been as a foundation of modern cosmology, 
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the triumph of GR is contingent upon one crucial caveat that the Universe is dominantly 
filled with an exotic energy component with negative pr essure (dubbed dark energy] which 



is be lieved to accelerate the Universe at present epoch (IRiess et al.l Il998l : iPerlmutter et al. 



19991 ). Although the cosmological constant A devised by Einstein himself is currently the 
paradigm of dark energy, the extreme fine-tuning of A required to explain the observations 
has provoked skepticism among the cosmologists about the validity of the standard ACDM 
cosmology (A+cold dark matter) based on GR. 

Recently, a fiurry of research has been conducted on modified gravity (MG) models, 
which attempt to modify GR on large scales so that the a pparent acceleration of the Uni- 
verse can be realized without introducing dark energy (see IClifton et al.ll2012l for a recent 
review). Among many different MG models that have been put forward so far, the f{R) 

l2ni_nh _ 



gravity (jSotiriou fc Faraonil l2010t Ide Felice fc Tsujikawa 



tentions after the recent stringent observational tests (IReyes et al.ll2010uWoitak et al.ll201ll ). 
Basically, the f{R) gravity generalizes the Ricci scalar R in the Einstein-Hilbert action for GR 
to a specific funct ion of R, which gives rise to an extra scalar degree of f reedom fji = df/ dR, 
dubbed scalaron (jSotiriou fc Faraonil I2OIOI : Ide Felice &: Tsujikawall2010l ). 



The scalaron produces an additional fifth force that modifies GR on large scales. In 
highly dense environments, howe ver, due to the chameleon rnechanism that effe ctively shields 
the fifth force, GR is recovered (iKhoury &: Weltmanll2004l : iLi &: Barrowll2007l ). For a given 
functional form of f{R), the formation of the cosmic structures and its clustering strength is 
determined by the magnitude of scalaron at the present epoch, |//jo|- The ACDM cosmology 
based on GR corresponds to the case of fm = 0, i.e., f{R) = A. The non-zero value of fm 
that is the essential feature of f{R) gravity, however, should be small enough to meet the 
observational constraints. 

Throughout this Pap er, we will focus ori the Hu-Sawicki f{R) gravity which has the 
following functional form (iHu fc Sawickill2007l ) 



fiR) 



-m 



ci{—R/m 



2\n 



C2{-R/m'^Y + r 



with two characteristic parameters of n and Ci/c2. Here, m = SnGpm/'^ with mean matter 
density pm at present epoch, and the free parameter Ci/c2 is related to the A density {fl\) 
and matter density [flm) parameters as Ci/c2 = GQ/i/Qm which mimics the evolution of the 
ACDM background. Regarding the other free para meter, n, we set its value at unity as in 
the previous works flOyaizull2008l : IZhao et al.ll201ll ). 



The observed abundance of galaxy clusters has constrained the present value of the 



scalaron in the Hu-Sawicki model to be IJrqI ^ 10 ^ (ISchmidt et al.ll2009l : iLombriser et al. 
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20101 ). while the solar system test p uts a more stringent constraint of \fm\ ^ 10^^ for the 
Milky way halo (iHu &: Sawickil 120071 ). When the scalaron meets this constraint, however, it 
is very difficult to distinguish the f{R) gravity model from the ACDM model based on GR 
since there is very little difference between the two models in their predictions for almost all 
observables. 

Here we suggest that it might be possible to discriminate f{R) gravity from GR by 
measuring the specific angular momentum distribution of the low-mass galactic halos located 
in low-density environments (i.e., the field regions) where the chameleon mechanism functions 
less effectively. We expect that due to the effect of the fifth force of the scalaron the low- 
mass field galactic halos in the f{R) gravity case would acquire higher specific angular 
momentum than those in the GR case. Whereas, the galactic halos located in high-density 
environments where the chameleon shielding mechanism effectively screens the fifth force 
would not s how any difference in the specific angular momentum distribution between the 
two models (jZhao et al.ll201ll : iLi et al.ll2012al ). 



Note that the expected enhancement in the specific angular momentum of the low-mass 
field galactic halos is essentially the same effect of the fifth for ce that e n hance s the velocity 
dispersion o f uns creened halos, which was first discussed by ISchmidtl ( 120101 ) and also by 
Zhao et al.l (1201 if ). The main difference is that for the case of the angular momentum what 
takes into account is only the velocity components perpendicular to the separation from the 
halo center on which the gravity act as a tidal torque forces. The enhanced gravity due to 
the presence of the fifth force leads to the enhanced tidal torque forces which in turn induces 
higher angular momentum of an unscreened halo. 

To investigate quantitatively how significantly the distribution of the specific angular 
momentum differs between the GR and the f{R) gravity cases, we use the high-resolution 
N-body simulations for the two models to determine the distributions of the specific angular 
momentum of the low-mass field galactic halos. Being conservative, we consider the case of 
|/i?o| = 10~^ as our fiducial f{R) gravity model, calling it "F6 model" from here on. The 
contents of this Paper are outlined as follows. In section [2] we describe the numerical data 
from the N-body simulations and present the results on the specific angular momentum and 
spin parameter distributions of the low-mass field galactic halos for the GR and the F6 cases. 
In section |3] we speculate a possibility of using the abundance of the low surface brightness 
galaxies as a probe of modified gravity. In section H] we discuss the results and draw a final 
conclusion. 
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2. TIDAL EFFECT OF THE F(R) GRAVITY 

2.1. Numerical data 

We hav e performed hig h-resolution A^-body simulations which implement the ECOS- 
MOG code (ILi et al.ll2012bl ) to compute the trajectories of 512^ dark matter particles in a 
periodic box of size 100 h~^Mpc on a side for the GR and F6 models. The key cosmolog- 
ical parameters that describe the initial conditions of each model are identically set at the 
WMAP7 values (n^ = 0.24, = 0.76, Vlb = 0,045, h = 0.73, as = 0.77, = 0.96) 
(IKomatsu et al.ll201ll ). The bound dark matter halos are identified from th e simu lation data 
by applying the AHF (Amiga's Halo Finder) code jKnoUmann fc Knebel I2OO9I ) which de- 
fines the halo mass M as the top-hat mass enclosed by the spherical radius R at which the 
overdensity reaches the threshold value of = 374 at z = 0. For the details of the simu- 



lations and the halo-identification procedures, see iKnoUmann fc Knebd (l2009l ) and iLi et al. 



Among the identified halos for each model, we select only those well resolved halos which 
consist of 100 or more dark matter particles {Np > 100). Then, we measure the physical 
angular momentum, J, of each selected halo as J = Yli ^ ^i^i ^ where mj, Vj and 
represent the mass, velocity and position of the i-th particle belonging to the halo. Here the 
positions of of the component particles, {ri}^^ , are all measured with respect to the center 
of mass of their host halo. We divide the selected halos into the field halos and the wall 
halos by applying the friends-of-friends (FoF) criterion which has been conventionally used 
in observational data analysis for the separation between the field and the wall galaxies. The 
usefulness of the FoF criterion lies in the fact that it requires no other information than the 
positions of the sample galaxies and thus th a t it ca n be readily appli e d to r eal observational 
data in practice. Very recently, IZhao et al.l (I2OIII ) and ICabre et al.l (120121 ) have developed 
more elaborate criteria with which the screened and unscreened halos can be separated. 
However, in the current work, we choose the simple FoF criterion for the practical purpose. 

We first measure the mean separation distance i among the selected halos, and then 
classify those isolated halos which do not have any neighboring halos within the threshold 
distance of hi as field halos, where h is the linkage length parameter. The rest of the selected 
halos are classified as wall halos which have at least one neighboring halo within residing 
in relatively higher density regions. The larger the value of h is, the more isolated the field 
galactic halos are. Here, we set the linkage length parameter at 6 = 0.2. We also restrict 
our analysis to the galactic halos which are less massive than a given mass threshold 
for which we consider two cases of = 10^^ and lO^^/i~^M0. In Tabled] we list the total 
number of the well- resolved halos with Np > 100 (Atotai), the mean halo separation i, the 
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numbers of the field galactic halos for two different mass thresholds. 



2.2. Distribution of the Specific Angular Momentum 

The specific angular momentum, j, of each field (wall) galactic halo is computed as j = 
J/M , in unit of [s~^km h~^Mpc]. Binning the magnitude of the specific angular momentum, 
j, and counting the number AA^^g, of the field (wall) galactic halos which belong to each 
j-bin, we determine the probability density distribution of the specific angular momentum 
of the field (wall) galactic halos as p{j) = ANg/{NgAj), where Ng is the total number of 
the field (wall) galactic halos and Aj represents the size of each j'-bin. 

Figure [1] shows the probability density distributions, p{j), of the field and wall galactic 
halos with the Jack-knife errors for two different cases of Mc in four separate panels. We 
construct eight Jack-knife resamples of equal size and then calculated aj as one standard 
deviation scatter of p{j) among the eight resamples. In each panel the solid and dashed lines 
represent the result for the GR and F6 cases, respectively. As can be seen, the galactic halos 
in the F6 model have higher specific angular momentum than those in GR. Furthermore, 
the field galactic halos with lower mass of M < 10^^ H^^Mq exhibit the largest difference 
between the two models. We interpret this result as follows. The enhanced gravity in the 
F6 model exerts enhanced tidal torque forces on the proto-galaxies in the field regions, and 
in consequence the field galactic halos for the F6 case to acquire higher specific angular 
momentum than for the GR case. 

To test the null hypothesis that there is no difference in the distribution of the specific 
angular momentum of the low-mass field galactic halos bet ween the GR and the F6 cases. 



we perform the two-sample KS (Kolmogorov-Smirnov) test ( iWall &: Jenkins! l2003l ). Figure 
[2] shows the cumulative fraction function of the specific angular momentum of the low-mass 
field galactic halos (M < 10^^ h'^M^) for the GR (solid line) and F6 (dashed line) cases. 
Measuring the KS statistic, i.e., the maximum distance between the two cumulative fraction 
functions is found to be 203, which rejects the null hypothesis at 99.9999% confidence level. 

Now that the difference in the distribution of j of the low-mass field galactic halos 
between the GR and F6 cases is found to be statistically significant, it is important to see 
whether the two cases have comparable mass distributions. For this, we first calculate the 
following three quantities: First, the probability density distribution of the logarithmic mass 
of the field galactic halos as p(lnM) = AA'g/(A'"gAlnM), whereA^g is the total number 
of the field galactic halos, AA"g is the number of the field galactic halos belonging to each 
logarithmic mass bin and A In M is the size of the logarithmic mass bin. Second, the number 
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densities of the field galactic halos per unit volume, dN/dlnM = AA'"g/(V^A In M) in unit of 
h^Mpc~^, where V represents the volume of the simulation box. Third, the mean mass of the 
field galactic halos belonging to each logarithmic mass bin, (M) = AM/ANg, where AM 
represents the sum of the masses of the field galactic halos belonging to each logarithmic 
mass bin. 

Figure [3] shows p{ln M) and dN / d In M in the bottom and top panels, respectively. In 
each panel the solid and dashed hues represent the GR and F6 cases, respectively. As can be 
seen, although the F6 model yields ~ 10% increment in the amplitude of the mass function 
of the field galactic halos, p(lnM) in the two models are almost identical to each other. 
Figure m shows the ratio of (M) for the F6 case to that for the GR case, which reveals that 
the difference in the mean mass between the two models is less than 1% in the whole mass 
range of the field galactic halos. The results shown in Figures [3}|1] prove that the two samples 
are comparable in the mass distribution and thus that the difference in j between the two 
models are not due to the mass-bias. 

To test if the difference in p{j) of the low- mass galactic halos between the GR and F6 
models shown in Figure [T] is really due to the effect of the fifth force but not due to any 
possible numerical flukes, we multiply the factor of (4/3)^/^ to the values of j of the low- 
mass field galactic halos for the GR case and determine its probability density distribution, 
p[(4/3)^/^j]. The gravitational tidal force in the low-density regions can be enhanced by a 
maximum factor of 4/3 for the F6 case. That is, the low-mass field galactic halos at their 
proto-galactic stages would experience enhanced gravitational tidal force in the F6 model. 
In consequence, the velocities that enter in the definition of the specific angular momentum 
would increase by a maximum factor of (4/3)^/^, which in turn leads to the enhanced specific 
angular momentum of the low-mass field galactic halos by the same factor. Therefore, it is 
expected that PGJ?[(4/3)^/^j] should match the probability density distribution of j of the 
low-mass field galactic halos for the F6 case, ppeU)- 

Figure [5] shows pGfi[(4/3)^''^j] (dotted line) and compares it with preU) (dashed line). 
As can be seen, the two distributions agree with each other very well, which confirms that the 
difference in p{j) of the low-mass field galactic halos is due to the effect of modified gravity 
which can be enhanced by a maximum factor of 4/ 3 in the low-den sity e nvironment s. This 



result is also consistent with the previous works of ISchmidtl (120101 ) and IZhao et al.l (120111 ) 
who showed that the effect of MG enhances the velocity dispersion of a unscreened halo by a 
factor of (4/3)^/^. The fact that the measurements of the velocity dispersion and the specific 
angular momentum yield an enhancement of the same factor for unscreened halos indicates 
that the velocity structure is unchanged by the modified forces at first order. 

Figure [6] shows the ratio of p{j) of the wall galactic halos to that of the field galactic 
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halos, Pwaii(j)/Pficid(j), for the GR and F6 cases. As can be seen, for GR, the ratio increases 
rapidly with j, while for the F6 case the rate of the increase of the ratio with j is much 
milder. The dominance of wall galactic halos in the high-j section is attributed to frequent 
merger events in high-density regions. As the merging rate is much lower in the field field 
regions for the GR case, it is natural to expect that most of the field galactic halos are biased 
toward the low j section. In contrast, for the F6 case, even though frequent merging events 
do not occur in the low-density regions, the unscreened fifth force has an effect of spinning 
up the field galactic halos, resulting in the increase of the relative abundance of the high-j 
field galactic halos and the decrease of Pwaii(j)/Pfieid(j) at the high-j end. 

One may also expect that the more field field galactic halos would acquire on average 
stronger effect of the enhanced gravitational tidal force since the degree of the chameleon 
effect is lower in the more field regions. To investigate how the probability density distribu- 
tion, p{j), depends on the degree of the isolation of the field galactic halos, we calculate p{j) 
repeatedly, varying the value of b from 0.1 to 0.3. Figure [7] shows the ratio of p{j) of the 
field galactic halos for the F6 case to that for the GR case, PfgU) / PgrU) , for three different 
cases of the linkage parameter b. As can be seen, the larger the value of b is, the larger the 
ratio, Pf6(j)/Pgr(j) is in the high-j section. This result is consistent with the picture that 
the chameleon screening of f{R) gravity becomes less effective in low-density regions. 



2.3. Distribution of the spin parameter 



The specific angular momentum of a galactic halo is commonly expressed in terms 
of the dimens ionless spin parame ter (jPeebled Il969l ). which is efficiently defined as A = 
?'/(2G'M-R)^'^^f Bullock et al.ll200ll . and references therein). N-body simulations found that 
the probability density distribution of this dir aensionless spin parameter is well approximate d 
by the following log-normal distribution (e.g. 

1 



p{X) 



Barnes fc Efstathiou 1987 : Bullock et al. I boOlh : 

In' (A/Ao 



exp 



2ax 



(2) 



Av27r(TA 

with Ao ~ 0.03 and ax ~ 0.5. N-body simulations also revealed that the shape of p(A) is 
almost independent of the initial condit ions of the Universe, always following the log-normal 
form with the same values of Aq and a\ ( iBarnes fc Efstathioulll987l : ISteinmetz fc Bartelmann 
19951 ). Furthermore, recent observations also showed that the spin parameters of the observed 



galaxies seem to be only very weakly correlated with the other physical propertie s of the 



1 . 

galaxies such as overdensity, mass, shape and internal structure (iDisney et al.l 120081 ). 



For each field galactic halo from our simulation, we compute the spin parameter and 
calculate the probability density distribution, p(A), for the case of the field galactic halos 
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with M < 10^^ H^^Mq. Fitting the numerically obtained spin parameter distribution, p(A), 
of the field galactic halos to Equation (|2]), we determine the best-fit values of Aq and ax 
for both of the GR and F6 cases. Figure [S] shows p{X) (square and circular dots) with the 
Jackknife errors and compares them with the best-fit log-normal distributions (solid and 
dashed lines) for the GR and F6 cases, respectively, in the top panel. As can be seen, p(A) 
is indeed well approximated by the log-normal distribution for both of the cases. Table [2] 
lists the best-fit values of Aq and ax for the two models. As can be read, for the GR case, 
the best-fit values of Aq ~ 0.03 and ax ~ 0.5 are consistent with the previous N-body results 



(IBuUock et al.ll200l[ ). while the F6 case yields larger values of Aq and ax as expected. 



Our result for the GR case confirms the claim of the previous works that the spin 
parameter distribution is insensitive to the halo mass scale, surrounding overdensity, and 
the background cosmology. Note that even though we have considered only the field galactic 
halos with mass less than 10^^ h~^MQ from a N-body simulation for a WMAP7 cosmology, 
we have obtained the lo g-normal distr i butio n of p(A) with Aq ~ 0.03 and ax ~ 0.5 which is 



identical to the result of iBuUock et al.l (120011 ) who used all dark matter halos from a N-body 



simulation with different cosmological parameters to derive p(A). 

Now that p{X) is found to differ significantly between the GR and F6 cases while it is 
insensitive to the background cosmology, the spin parameter distribution of the low-mass 
field galactic halos should be a sensitive indicator of underlying gravity. The bottom panel 
of Fig. [H] shows the ratio of the probability density distribution of A for the F6 case to that 
for GR, Pf6(A)/pgr(A). As can be seen, the ratio exceeds unity at A ~ 0.05 and increases 
sharply in the high-A section (A > 0.05). The cumulative probability of P(A > 0.05) are 
found to be 0.22 and 0.31 for the GR and F6 cases, respectively. 



IMPLICATION ON THE LOW-SURFACE BRIGHTNESS GALAXIES 



The results presented in section 12.31 has an interesting possible implication on the abun- 
dance of the low surface brightness galaxies (LSBGs) whose apparent fiuxes are an order of 
magnitude lower than the night sky. The LSBGs have b een found to occupy a significant 
fraction (~ 30 — 50%) of the total field galaxy po pulation (llmpey fc Bothunlll997l . and refer- 
ences therein). Very recently, iGeller et al.l (120121 ) determined the galaxy luminosity function, 
using a magnitude limited sample of the galaxies with i?-band magnitude less than 20.6 from 
the Smithsonian Hectospec Lensing Survey (SHELS) which is a dense redshift survey cov- 
ering a four degree square region at redshifts 0.02 < z < 0.1. They found that the slope of 
the luminosity function in the faint-end increases up to 15.2 from 13.2 when the LSB drawf 
galaxies in the field are included. 
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Regarding the origin of the LSBGs, several theoretical studies suggested a hypothesis 
that the only difference of the LSBGs from t he ordinary galaxies is their biased formation in 
the fast-spinning dark halos with A > 0.05 ( llmpey fc BothunI 119971 . and references therein). 
The faster spinning motion of a host halo results in the decrements of the surface density of 
its disk gal axy, which in turn is directly proportional to the surface brightness. To test this 
hypothesis, Ijimenez et al.l ( 1l998l ) constructed a theoretical model under the assumption that 
the dark halos where the LSBGs reside have relatively higher specific angular momentum 
with the spin parameters larger than 0.05, and compared the predictions of their model with 
the observed properties of the LSBGs. According to their results, the theoretically predicted 
properties of the LSBGs such as colors, metalli ticty and spect r a are all in excellent agree- 
ments with the observe d values. The results of Jimenez et al.l (1998 ) have been confirmed 
by the follow-up works (jBoissier fc Prantzos 200ol . 



2001 



Boissier et al.l 120031 ) 



Now that the LSBGs are found to occupy a significant fraction of the field galaxy 
population (over 30%) and located in the fast-spinning dark halos, a fundamental question is 
whether or not the standard ACDM cosmology can produce such abundant fast-spinning dark 
halos in the field regions. At face values, our results then imply that MG, if existent, would 
produce more LSBGs in the field regions. The difference between the relative abundances of 
the LSBGs in the GR and F6 cases (22% and 31%, respectively) is definitely significant even 
though the F6 model is very similar to GR in its predictions for most other observables. 

Besides, it has been known that the spin parameters of dark matt er halos depend only 
very w eakly on the initial conditions of the Universe. For instance, iBarnes fc Efstathiou 
(119871 ) have shown by N-body simulation s that the spin parameters of d ark halos are un- 
correlated with the initial overdensities. ISteinmetz fc Bartelmanru ( 119951 ) have also found 
through N-body simulations that there is very weak correlation between the spin parameter 
A and the initial rms density fluctuation a, concluding that the s pin paramete r s shou ld be 
independent of the type of the initial power spectrum. Recently, iDisney et al.l (120081 ) have 
found by analyzing the observational data that the spin parameters of the local disk galaxies 
are uncorrelated with the other physical properties of them. 

Given our result that the F6 model yields more dark halos with A > 0.05 and the 
previous results that the spin parameter is almost independent of the initial conditions of 
the Universe, the relative fraction of the LSBGs that are believed to reside in halos with 
A > 0.05 could be an optimal indicator of the effect of MG. 
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4. SUMMARY AND CONCLUSION 

Analyzing the numerical data from the high-resolution N-body simulations for the Hu- 
Sawicki f{R) gravity with \fRo\ = 10~® (F6 model), we have shown that the probability 
density distribution, p{j), of the specific angular momentum of the low-mass field galactic 
halos with M < 10^^ h'^M^ in the F6 model is significantly shifted toward the higher j 
section than in the GR case. The null hypothesis that there is no difference in p{j) between 
the GR and the F6 models is rejected via the KS test at the 99.999% confidence level. 

We have also determined from the numerical data the probability density distributions 
of the dimensionless spin parameter, p{X), for the GR and the F6 models and found that the 
high-spin field galactic halos with A > 0.05 occupy 21% and 31% of the total field galactic 
halo populations for the GR and the F6 models, respectively. Given this result and under 
the assumption that the host halos of the low-surface brightness galaxies (LSBGs) have spin 
parameters larger than 0.05, we speculate a possibility of discriminating the f{R) gravity 
from the GR with the relative abundance of the LSBGs in the field regions. 

In order to use the LSBGs as a test of gravity, however, it will be first necessary to as- 
sess the p ossible systematic involved in modeling the surface brightness of the disk galaxies. 



Recently, iDavis et al.l ( 12012| ) recently claimed that the dwarf void galaxies could be signifi- 
cantly brighter in MG models with chameleon mechanism since the unscreened MG in void 
regions would play a role of enhancing the intrinsic luminosity of the main sequence stars. 
The situation could be subtle when other factors are taken into account. For example, the 
unscreened MG would not only brighten up the stars, but also change the dynamical prop- 
erties of the host halos. These effects need to be taken into account to relate our predictions 
for dark matter halos to the abundance of the LSBGs. We plan to work on this issue and 
hope to report the result elsewhere in the near future. 
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Fig. 1. — Probability density distributions of the specific angular momentum of the field 
and wall galactic halos in the left and right panels, respectively. The top and bottom panels 
correspond to the case that the galactic halos with masses of M < 10^"^ H^^Mq and of 
M < 10^^/i"^Mq, respectively. In each panel, the solid and the dashed lines correspond 
to the standard GR+ACDM and F6 case, respectively. The errors are obtained using eight 
Jackknife resamples. 
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Fig. 2. — K-S test cumulative fractional plot of the specific angular momentum of the field 
galactic halos for the two models. 
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Fig. 3. — Number and probability densities of the galactic halos as a function of their 
logarithmic mass in the top and bottom panel, respectively. 
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Fig. 4. — Ratio of the mean mass of the galactic halos in the F6 model to that in the GR 
model. 
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Fig. 5. — Same as the top-left panel of Figure [H but plotting also the converted probability 
density distribution of j (dotted line) obtained by multiplying the values of j for the GR 
case by a factor of (4/3)"'^/^. 
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Fig. 6. — Ratio of the probability density of the field galactic halos to that of the wall galactic 
halos for the two models. 
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Fig. 7. — Ratio of the probability density of the field galactic halos in the GR model to that 
in the F6 model for the three different cases of the linkage length parameter. 
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Fig. 8. — (Top panel): Probability density distributions of the spin parameter of the field 
galactic halos for the GR and F6 case as square and triangle dots, respectively. The sohd 
and dashed lines correspond to the best-fit log-normal distributions for the GR and F6 case, 
respectively. (Bottom panel:): Ratio of p(A) for the F6 case to that for the GR case. The 
dotted line corresponds to the fiat case. 
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Table 1. model, # of the well-resolved halos with > 100 (iVtotai); the mean halo 
separation, # of the field galactic halos for the two different mass thresholds 



model 




e 


-^field.I 


-^field.II 








(M < 10^2 h-^M^) 


(M < IO^i/i-^Mq) 


GR 


43396 


2.846 


24114 


11252 


F6 


47517 


2.761 


26297 


12299 
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Table 2. model, two best-fit values of the spin parameter distribution 



model 


Ao 




GR 


0.033 


0.556 


F6 


0.038 


0.569 



